BI-ORTHOGONAL POLYNOMIALS ON THE UNIT CIRCLE, 
REGULAR SEMI-CLASSICAL WEIGHTS AND INTEGRABLE 

SYSTEMS 

P.J. FORRESTER AND N.S. WITTE 

Abstract. The theory of bi-orthogonal polynomials on the unit circle is de- 
veloped for a general class of weights leading to systems of recurrence rela- 
tions and derivatives of the polynomials and their associated functions, and to 
functional-difference equations of certain coefficient functions appearing in the 
theory. A natural formulation of the Riemann-Hilbert problem is presented 
which has as its solution the above system of bi-orthogonal polynomials and 
associated functions. In particular for the case of regular semi-classical weights 
on the unit circle w(z) = IljLiC- 2 — 2 j(*)) Pj > consisting of m £ Z>o finite sin- 
gularities, difference equations with respect to the bi-orthogonal polynomial 
degree n (Laguerre- Freud equations or discrete analogs of the Schlcsinger equa- 
tions) and differential equations with respect to the deformation variables Zj (t) 
(Schlesinger equations) are derived completely characterising the system. 



1. Introduction 

The unitary group U (N) with Haar (uniform) measure has eigenvalue probabil- 
ity density function (see e.g. [23 Chapter 2]) 

(L1) (2nWM II l**-*il 2 ' *«:=e*'eT, (9; G {—%, tt], 

where T = {zgC:|z| = 1}. Our interest is in averages over U G U(N) of class 
functions w(U) (i.e. symmetric functions of the eigenvalues of U only) which have 
the factorization property YliLi w i z i) f° r { z i? ■ • ■ ? z n} G Spec(J7). Introducing the 
Fourier components {w;}z e z of the weight w(z) — J>2i^_ 00 w i zl > due to the well 
known identity |45) 

N 

(1.2) (l[ w ^) u(N) = det K-i]i,j=i,-,JV, 

i=i ( ' 

we are equivalently studying Toeplitz determinants. However we are interested in 
the situation where the weights are not necessarily positive or even real valued 
w(z) ^ w(z) for z£T, where the bar denotes the complex conjugate, and conse- 
quently the Toeplitz matrices are non-Hermitian w n ^ w- n - The motivations for 
studying these types of weights comes from many diverse applications, in particular 
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the gap probabilities and characteristic polynomial averages in random matrix the- 
ory, the spin-spin correlations of the square lattice Ising model, the density matrix 
of one-dimensional systems of impenetrable bosons and probability distributions for 
various classes of non-intersecting lattice path problems. All the above applications 
correspond to special cases of the weight 



where fi, w = U\ +itU2 are complex parameters and £, t = e l ® are complex variables. 
We have previously shown that l|1.2f) with weight l|1.3f) can, as a function of t, 
be characterised as a r-function for the Painleve VI system, and thus be expressed 
in terms of a solution of the sixth Painleve equation. Here we will develop theory 
which allows for a characterisation of a general class of Toeplitz determinants as 
a function of the discrete variable N. In the particular case of the weight (|1.3|) 
the resulting difference equations in N will be studied in a companion paper |25| 
where it is shown they are equivalent to the discrete fifth Painleve equation dPy 
associated with the degeneration of the rational surface D^ 1 — > . 

To characterise the unitary group averages l|1.2l) we have found it necessary to 
substantially develop the theory of bi-orthogonal polynomial systems on the unit 
circle {(/> n (z), </>n(z)}^0' generalising the Szego and Geronimus theory applying to 
orthogonal polynomial systems. They arc defined by 



Issues relating to the existence of such systems and their basic properties is taken up 
in Section 2. A study of bi-orthogonal polynomial systems on the unit circle was 
begun by Baxter [111 112) where their elementary properties were elucidated and 
the connection between the absolute convergence of the trigonometric expansion 
of \ogw(e l6 ) and absolute convergence of certain coefficients of the bi-orthogonal 
polynomials was investigated. Sometime later other useful elements of the theory 
of orthogonal polynomial systems were added by Jones, Njastad and Thron in their 
study j35| of the trigonometric moment problem and related quadrature formulae 
on the unit circle, namely the associated polynomials and the Hilbert transform 
of the weight. The task of determining the functional, difference and differential 
equations for orthogonal polynomial systems on the unit circle was initiated in 
|33| but was not completed nor extended to cover the bi-orthogonal situation. In 
extending the theory in this direction we have sought to cast it in a way so that the 
appearance is as close to the Szego- Geronimus theory as possible, and indeed we 
find that virtually all formulae from the older theory can be taken over but now the 
complex conjugated variables have to be re-interpreted as independent variables. 
We wish to emphasise that the works described above, which are directly relevant 



(1.3) 




1, 9 S (-7T,7T - (j)) 

€ (n - <j>,ir] 



(1.4) 
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to the present study of certain classes of weights, constitute only a small fraction 
of existing literature on orthogonal polynomial systems on the unit circle. The 
themes taken up in this broader body of work, for example as recounted in |44) . are 
really rather different problems and are formulated for the most general types of 
real, positive measures on the unit circle. Finally we note that to a large extent the 
task of characterising the analogous Hankel determinants has been completed for 
weights defined on the real line and this could be founded upon the standard theory 
of orthogonal polynomial systems. In particular we note the works of Bauldry 10 , 
Bonan and Clark ^B], Belmehdi and Ronveaux and Magnus [3SI El EM] where 
in the last three works the differential and difference structures have been clearly 
revealed and linked to isomonodromy preserving deformations, culminating in the 
identification of the simplest nontrivial case with the sixth Painleve system. 

In our formulation of the problem we start with a weight, possessing certain 
features to be made precise later, which we take as given and wish to calculate the 
Toeplitz determinants via auxiliary quantities arising in the bi-orthogonal polyno- 
mial theory. To this end we have derived closed systems of differential relations for 
the polynomials {4>n(z)}^ =0 , their reciprocal polynomials {^(^Oj-^Loi an< ^ associ- 
ated functions {e n (z)}*^L , {^n( z )}ri=o m Proposition ^. 61 In the notation of i 
and Corollary 2.3 let 



Y n (z;t) 

We show 




e n (z)/w(z) 

-£*n(z)/w(z) 



—Y - AY 

dz 



where the entries in the matrix A n are parameterised by four coefficient functions 
Q n {z), 0*(z) , 6„(z) , 6* (z) in (|2.71|) . and complete sets of difference and func- 
tional relations for these coefficient functions are given in Proposition ^ . 71 and Corol- 
larv !2.2l We also formulate a 2 x 2 matrix Riemann-Hilbert problem in Proposition 
12. 81 for general classes of weights which parallels the case for orthogonal polynomials 
on the line |^ |^ 1221 IT%| and whose solution is simply related to Y n . 

For our particular applications the weights are members of the regular semi- 
classical class 

m 

w{z) = \{{z-z 3 {t))P^ Pj ec, 

3=1 

with an arbitrary number m of isolated finite singularities located at zj(t). These 
are also known as generalised Jacobi weights. Using the projection of the unit circle 
onto the interval [—1,1] a system of orthogonal polynomials on the unit circle is 
equivalent to two related systems of orthogonal polynomials on this interval. The 
generalised Jacobi orthogonal polynomial systems with support [— 1 , 1] have been 
studied from a number of points of view. For example the Stieltjes electrostatic 
problem was generalised to include an arbitrary number of fixed charges in |27j . 
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Uniform asymptotics of the orthogonal polynomials were first derived in [5], [Hj and 
this was extended to their derivatives in 0], 0H|> |3Z| an d |49| . Further properties 
of the polynomials, such as their upper bounds, were investigated in |19| . |43) and 
estimates of the associated functions and polynomials were made in |17| . Interest in 
issues relating to quadrature problems were taken up in |2()| where lower and upper 
bounds of the Christoffel function were found and quadrature inequalities were 
derived in |42|. The zeros of the orthogonal polynomials have also been studied, in 
particular the spacing of consecutive zeros in |20] and their asymptotic formulae in 
|46j . We wish to mention too a work similar in spirit to our own |26| . where the 
orthogonal polynomial system {p n (t)}^Lo was defined by 

/ dtw(t)p n (t)t k = 0, fc = 0,...,n-l. 
Jc 

In this work the weight had m — 3 distinct finite singularities of the regular semi- 
classical class and its support was a closed curve C enclosing all the singularities. 
Recurrence relations in n = deg(p„) for the three-term recurrence coefficients were 
found and heuristic arguments were given for the asymptotic behaviour of these as 
well as the polynomials themselves. 

A key feature of regular semi-classical weights is that 

1 d _ 2V(z) 
w{z) dz { ' W(z) ' 

where V and W are polynomials such that degV(z) < m,degW(z) = m. The 
coefficient functions for regular semi-classical weights are polynomials of z with 
bounded degree, degf2„(z) = degfi* (z) = m — 1, deg6„(z) = degO* (z) = m — 2 
(see Proposition I3.I|I . In addition evaluations of these functions at the singular 
points satisfy bilinear relations (see Proposition 13. 2fl which lead directly to one of 
the pair of coupled discrete Painleve equations. Deformation derivatives of the 
linear system of differential equations above with respect to arbitrary trajectories 
of the finite singularities are given in Proposition 13 . 31 which can summarised as 

d Yn . = Bn Y n = J 5oo -V^l Zj .iy n; where A n = V 

dt | U z - Zsdt J z ~ Z] 

and consequently systems of Schlesinger equations for the elements of A n j (or the 
coefficient functions evaluated at Zj) are given in (|3.77I3.79|I . It is quite natural that 
systems governed by regular semi-classical weights preserve the monodromy data of 
the solutions Y n about each singularity Zj with respect to arbitrary deformations. 

We mention at this point that the full definition of a regular semi-classical weight 
given below in Defmition l3.Il is restrictive, and has been generalised by relaxing some 
of the conditions in a series of works , jSH] > . In these works the orthogonal 
polynomial systems were characterised by integral representations of semi-classical 
linear functionals with respect to certain paths in the complex plane. An irregular 
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semi-classical weight arising under less restrictive conditions can be recovered from 
a particular regular semi-classical weight through limiting processes involving the 
coalescence of singular points in the same way as the fifth to the first Painleve 
systems are recovered from the sixth. The differential and difference structures of 
orthogonal and bi-orthogonal polynomial systems with polynomial log-derivatives 
of the weight (deg(W) = 0, deg(V) > 2) defined on M have been studied ^Uj Q5| in 
the context of matrix models. Analogous structures for a bi-orthogonal polynomial 
system on the unit circle with a simple Laurent polynomial log-derivative for the 
weight were found in [31], motivated by applications to unitary matrix models. 

In Section 2 we derive systems of differential-difference and functional relations 
for bi-orthogonal polynomials and associated functions on the unit circle for a gen- 
eral class of weights and formulate the Riemann-Hilbert problem. In Section 3 we 
specialise to regular semi-classical weights and derive bilinear difference equations. 
In addition we calculate the deformation derivatives of the bi-orthogonal polyno- 
mial system, derive a system of Schlesinger equations and show the deformations 
are of the isomonodromic type. 



We consider a complex function for our weight w(z), analytic in the cut complex 
z-plane and which possesses a Fourier expansion 



where zeDcC and T denotes the unit circle |£| = 1 with £ = e l9 ,8 G (— 7r, tt]. 
Hereafter we will assume that z^w(z), z^w'(z) 6 L 1 (T) for all j E Z. We will also 
assume that the trigonometric sum converges in an annulus D = {z 6 C : Ai < 
\z\ < A2} and TcO. Until we arrive at our specific weights of interest in Section 
3, namely the regular semi-classical class, we will formally assume convergence 
holds. The doubly infinite sequence {wfe}^L_ DO are the trigonometric moments of 
the distribution w(e t6 )d9 /2ir and define the trigonometric moment problem. Define 
the Toeplitz determinants 



2. Bl-ORTHOGONAL POLYNOMIALS ON THE UNIT CIRCLE AND 



Riemann-Hilbert Problem 




00 



(2.1) 




det {w- e+j -k] 



0<j,k<n—l ' 



(2.2) 




where e will take the integer values 0, ±1. 
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We define a system of bi-orthogonal polynomials {4> n (z), <j) n (z)}^' =0 with respect 
to the weight w(z) on the unit circle by the orthogonality relation 

(2.3) f -^-^(0^(0^(0 = W 

This system is taken to be orthonormal and the coefficients in a monomial basis 
arc defined by 

n 

(2.4) <j> n {z) = K n z n + l n z n ~ x + m n z n - 2 + . . . + <f> n (0) ^J^c^z 3 , 

3=0 

n 

(2.5) 4> n {z) = R n z n + l n z n - x + rh n z n - 2 + ... + 0„(O) = ^ c n ^z\ 

3=0 



where R n is chosen to be equal to n n without loss of generality (this has the effect of 
rendering many results formally identical to the pre-existing theory of orthogonal 
polynomials). Notwithstanding the notation c n j in general is not equal to the 
complex conjugate of c n j and is independent of it. We also define the reciprocal 
polynomial by 



n 

(2.6) <&{z):=z n t n {l/z) = Y,*njz 

3=0 



The bi-orthogonal polynomials can be defined up to an overall factor by the or- 
thogonality with respect to the monomials 

(2-7) / 2Sc w(C)0Tl(C) ^ = ° °^3<n-l, 

whereas their reciprocal polynomials can be similarly defined by 
(2.8) ^_g_ w(C ) C(C xJ = l<j<n. 

The linear system of equations for the coefficients c n j , c n j arising from 



(2-9) c n , n J^^w(C)MCK 



dC - 0< to < n - 1 

1 to = n 



dC - - 00<m<n-l 

(2-10) c„,„ _^ w (C)r^(C) ={ ~ ~ 

2m( h m = n 
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has the solution 



det 



(2.11) 



Wo 

W n -l 
V W n 



W-l 
Wo J 



det 



w 



\w n 



W-n\ 

w J 



det 



(2.12) 



w 

W-n+1 
\ W- n 



Wl 
W J 



( 



det 



w 



\W-r 



Wn 



w J 



and in particular one has the following results. 

Proposition 2.1 (|llj). The leading and trailing coefficients of the polynomials 
4>n{z), 4> n {z) are 

1 1° 

(2.13) 
(2.14) 



CnO = 0n(O) - (-1)" — C n o = MO) = (-1)" 1 '•• 1 



, TO ■ 
v « 1 n+l 



K n I„+l 



Proposition 2.2 (11 ). The bi-orthogonal system {0„(z), </>* (z)}^L exists if and 
only if I® 7^ for all n > 1 . 

Remark 2.1. We shall see that a failure of this condition can and generically must 
occur in the case of regular semi-classical weights which contain deformation pa- 
rameters z\ , . . . , z m and this is precisely the condition for a movable singularity (in 
this case a pole) in the dynamics with respect to the Zj. 

A consequence of the above solutions are the following detcrminantal and inte- 
gral representations for the polynomials, 

f w a ... W- j . . . W- n \ 



(2.15) 



(2.16) 



-l) n K, 



W n -1 

V i ... 

JXCXC-*)] 



Wn-j-X 
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(2.17) 



(2.18) 



( z ) = -fi det 



W-n+1 Z n \ 



W-J + l z J 



\ W n ... Wi 

/°mc)(i-0] 



ij 



The system is alternatively defined by the sequence of ratios r n — 0„(O)/k„, 
known as reflection coefficients because of their role in the scattering theory formu- 
lation of OPS on the unit circle, together with a companion quantity f n — </>„(0)//c ra . 
As in the Szego theory 0^1 r„ and f n are related to the above Toeplitz determinants 

by 



(2.19) 



rn = (-1)' 



B\w\ 



= (-1)'' 



Jn>] 



The Toeplitz determinants of central interest can then be recovered through the 
following result. 

Proposition 2.3 ( 12 ). With the convention I§ = 1 the sequence of {In}^Lo 
satisfy 



(2.20) 



= 1 - r n r n , n>l. 



(» 2 

subject to the condition r n f n 1 for all n > 1. 

Fundamental consequences of the orthogonality condition are the following cou- 
pled linear recurrence relations. 

Proposition 2.4 (QUE]). 



(2.21) 
(2.22) 



K n (j) n+1 (z) = K n+1 Z(j) n (z) + <fi n+1 (0)(j>n(z), 
K n (j)* n+1 (z) = K„+10,*(z) + 0„+i(O)z0„(z). 



One finds three-term or second order recurrences for the uncoupled recurrence 
relations 

(2.23) 

(0)</Wi(z) + K n _i0 n+1 (O)^0 n _i(z) = [k„0„+i(O) + (O)z]0„(z), 

(2.24) 

K n cj) n {0)<j)* n+l (z) + k„_i0„ + i(O)z0,*_ 1 (z) = [K n cj) n+ i(0)z + K„ + i0 n (O)]0* (z). 

The analogue of the Christoffel-Darboux summation formula is given by the 
following result. 
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Proposition 2.5 ( 12 ). 

(2.25) 2^,9j{ z m\.Q = T^~7 : 

3=0 Z ^ 



^ = 0n+l(^X+l(C) ~ </>n+l(z)</>n+l(C) 

1 - z( 

for zQ ^ 1. -ffere 

(2-27) 0* (C) = C"^(VC)- 

Equations (|2.23I2.24|) being second order linear difference equations admit other 
linearly independent solutions ip n (z), ip*(z), and we define two such polynomial 
solutions, the polynomials of the second kind or associated polynomials 

(2.28) ip n (z) := f t^7t-^w(C)[0„(C) - <t> n (z)], n > 1, Vo : = «oW = 1/kq, 
and its reciprocal polynomial ^>*(z). The integral formula for -0* is 

(2.29) r n (z) ■=- I 9 dC . C ? + Z w(C)[z n MC) - <t>* n (z)l n>l, ^ := l/« . 

A central object in the theory is the Caratheodory function, or generating function 
of the Toeplitz elements 

(2-30) F(z) := / A,i±£u,(C), 

which has the expansions inside and outside the unit circle 

wq + 2^,^=1 w kZ k , if |z| < A min < 1, 

-wo - 2 X)fcli w -kz~ k , if |z| > A TOax > 1. 



(2.31) F(z) 



Having these definitions one requires two non-polynomial solutions e n (z),e*(z) to 
the recurrences and these are constructed as linear combinations of the polynomial 
solutions according to 

(2.32) e n (z) := ^ n (z) + F(z)q> n (z) = / -^-jj + Z w(Q<p n (0, 



T 27TZC C - * 

e* n ( Z ) ■= r n (z) - F(zmz) = -z» / dC C+Z w(0M0, 

(2-33) =- - / ^^"(CMUO- 

These have integral representations analogous to (|2.16I2.18|I 



(2.34) -£e n (z) = z 



2 nv 7 CiKO] ' 

(2.35) ^;(^(-r J ° +lWC)(c " r11 



2 nw v ; ^+iK0] 
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for \z\ ^ 1, which are particular cases of more general moments of the characteristic 
polynomial considered by Ismail and Riidemann |32| . 

Theorem 2.1 ( |2H1 ; EHj i EOj i EH] ) • VViW: satisfy the three-term recurrence 

relations \2.2'A \2.2$ and along with e n (z),e*(z) satisfy a variant of \2. 2112.22]) 

namely 

(2.36) K n e n+1 (z) = K n+ \ze n {z) - 0„+i(O)e* (z), 

(2.37) K„e* +1 (z) = K„+ie* (z) - 4> n+ i(Q)ze n (z). 

Theorem 2.2 f|28|). The Casoratians of the polynomial solutions (j) n , <f>* n , ip n , ip„ 
are 



(2.38) 
(2.39) 
(2.40) 



(f) n+ i(z)ip n (z) - 4> n+ i(z)(j) n (z) = 4> n+ i{z)e n {z) - e n+ i(z)tp n (z) = 2 ^ n+1 ^ z r 



& n (z)Tp*(z) + Mz)K( z ) = <Pn(zK(z) + en{z)<l>*(z) = 2z n . 

Further identities from the Szego theory that generalise are those that relate 
the leading coefficients back to the reflection coefficients 

(2.41) 4 = ^+^(0)^(0), 

(2.42) hL = h±L + rnfn _ u 

K n K n— 1 



(2.43) — = 

^n ^n — 1 



, _ ln-2 

r n -2 + r n ~i 

Kn-2 



Some useful relations for the leading coefficients of the product of a monomial and 
an bi-orthogonal polynomial or its derivative are 

(2.44) 

Z(j) n {z) = -—^—4) n+ l(z) + (— - hl±L \ (j) n ( Z ) 
Kn+l \ K n K n+1 / 

ln+1 In \ m n m-n+1 K r< 



K n— 1 \ K n+l K n / K n — 1 K n+1 K n — 1 



4>n-l(z) + 7T n _ 2 , 



Z 2 (/) n (z) = —^—(t> n +2{z) + ( — — j n+1 ( z ) 

Kn+2 \Kn+l Kn+2 K n +i / 

\ ln+1 f ln+2 In \ . Tn n TU n+ 2 1 , / \ , 

+ < H > <Pn(z) + 7T„_1, 

L tn+1 \K n+ 2 Kn J K n K n+2 J 
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<t>' n (z) = n — <f> n -i(z) + 7T„_ 2 , 

K n~l 

(2.45) z <t> n ( z ) = n( t>n{z) — 4>n-x{z) + 7T n _2, 

z 2 4>' n (z) = n-^—(j) n+ i{z) + |(n - 1)— - n^-^- \ (j) n {z) + n n -u 

where ' denotes the derivative with respect to z and where 7r„ denotes an arbi- 
trary polynomial of the linear space of polynomials with degree at most n. These 
identities can be verified directly. 

We will require the leading order terms in expansions of 4> n {z), 4>n( z )> e n{z), e* (z) 
both inside and outside the unit circle. 



Corollary 2.1. The bi- orthogonal polynomials 4> n {z) , <t>ni z ) have the following ex- 
pansions 
(2.46) 

I (o) + — ( (O) + 0„(O)r„_!)z + O(z 2 ) |z|<l, 

<t>n{z) = < 

{n n z n + l n z n - x + 0(z"- 2 ) |z| > 1, 

(2.47) 

|K n + r„z + 0(z 2 ) \z\ <1, 

" (Z) ~ I 0™(O)z" + _^( Kn ^„_ 1 (0) + ^.(OjZn-iJz"- 1 + 0(z"~ 2 ) |z| > 1, 

k K n -1 

whilst the associated functions have the expansions 
(2.48) 

' ' ' 0(z n + 2 ) 



l -e n (z) 



0n+l(O) _ x 

z 



|z| <1, 



K 2 0„ +2 (O) </>„+! (0) I 



K n+1 K n+2 



K n+1 ^n+1 



n+1 \ z - 2 



+0(z- 3 ) 



|z| > 1, 



(2.49) 



0n+l(O) _„ + i 
Kn+1 

1 Z 



K 2 </>„ +2 (0) </>„+! (0) I 



K n+1 K n+2 



"+ 1 ) z n+2 

+0(z™+ 3 ) \z\ < 1, 



z- 2 + 0(z- 3 ) |z|>l, 



«n+l \ K «+2 K n+l K n+2 

Proof. The second line of Ij2.46|l and the first line of (|2.47|) follow from the defini- 
tions. For the remaining lines of these two formulae we use 

(2.50) Kn-l^n(O) (0) + </>„(0)Z 

n — 1 j 

which results from differentiating ({2.21(1 and setting z = 0. The first line of (|2.48(1 
and the second line of (|2.49|l can be derived by employing the uncoupled recurrence 
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relations 1(2.23(1 and ((2.24(1 respectively. The remaining lines of these two formulae 
can be found by using the identity 

(2.51) z<t> n {z) = J^cj> n+1 {z) -^r^(0)^(4 

which in turn follows from combining 1)2.21(1 and ((2.22(1 . □ 

The z-derivatives or spectral derivatives of the bi-orthogonal polynomials in 
general are related to two consecutive polynomials [33| and we extend this to all 
cj> n (z), </>Jj(.z), e n( z )i e n( z ) with the following parameterisation. 

Proposition 2.6. The derivatives of the bi-orthogonal polynomials and associated 
functions are expressible as linear combinations in a related way (' := d/dz), 

(2.52) W{z)4>' n {z) = Q n {z)4> n+ i(z) - (n n (z) + V{z))<f> n {z), 

(2.53) W(z)^(z) = -e*„(z)^ +1 (z) + (n*(*) - V(z))<t>* n (z), 

(2.54) W(z)e' n (z) = Q n {z)e n+l {z) - (0„(z) - V{z))e n {z), 

(2.55) W{z)e*;{z) = -Q* n {z)e* n+1 {z) + (n* n (z) + V(z))e* n (z), 
with coefficient functions W(z),V(z) independent of n. 

Proof. The first, 1(2.52(1 . was found in (33j where the coefficients were taken to be 
(their notation A n , B n should not be confused with our use of it subsequently) 
k„_i0„ + i(O) zQ n (z) 



(2.56) A n = - 

1 / r A.. , , (n\ , , i 

e„{z) 



The other differential relations can be found in an analogous manner. □ 



(2.57) Bn = — \n n (z)+V(z) 



-z 



The coefficient functions 0„(z), il n (z), 9* (z), f2* (z) satisfy coupled linear re- 
currence relations themselves, one of which was reported in |33j . The full set are 
given in the following proposition. 

Proposition 2.7. The coefficient functions satisfy the coupled linear recurrence 
relations 

(2.58) U n (z) + U n -i(z) - H z <dn{z) + (n - 1) = 0, 

V (p n (0) K n ) z 

V q>n{0) K n J 

K n <f> n+2 (0) K„_i(/) n+ i(0) 4> n+ i(0)W(z) 

K n+ 10„ + l(O) Kn0n(O) <Pn(0) Z 

(2.60) ft„(z) + n n _ 1 (2) - + - z 6„(z) - n = 0, 
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(2.61) 



Vl+l . Yn+1 



(0) 



K n <p n (0) 

«rt^n+2(0) 
K„+10„+l(O) 



K n (f> n {0) K„ Z 



= 0, 



(2.62) 



^n+2(0) K n+2 \ , , K n+ i . . 

77^ + 2 B„+i(z) + (z6 n (z) - 6„(z)) = 0, 



(2.63) n n ( z ) - n n+1 (z) + 



Kn+2 { . 0n+l(O) 0^+2 (0) \ , 

2 H I On+i^Zj 



«7»+l 



tn+1 Kn+2 



(2.64) n; l+1 (z) + n n (z) 



K n +lK n K n Z 

K n +2 0n+2(O) 
Kn+1 0n+l(O) 

K n +1 



* ©;+i(*) 



(*e n (*) - e;(z)) - ^ = o, 



(2.65) n„(0) - o n+l (2) + — i + z e n+l (z) 

K n +l V. Kn+l Kn+2 / 

(/)„+! (O)0„+i(O) Kn+l 
H zSn(z) ©n( z ) = 0- 

Proof. The first (|2.58|) was found in by a direct evaluation of the left-hand 
side using integral definitions of the coefficient functions, however all of the re- 
lations follow from the compatibility of the differential relations and the recur- 
rence relations. Thus H2.58I2.59H follow from the compatibility of l|2.52|l and (|2.23f> . 
$n5ty . $nn\ follow from (|i2~53|) and (E2H>, follow from the combination 

of (|152li>.5^j> and ffiZIJl . and (|2\ 6411651) follow from the combination of (|152ll5:]|) 
and QTZfy . D 

Remark 2.2. The relations given above are obviously not all independent, as for 
example we note that (|2.58|l can derived from (|2.62l) with the use of (I2.67f) below. 

Corollary 2.2. Some additional identities satisfied by the coefficient functions are 
the following 



(2.66) 
(2.67) 
(2.68) 



<M+1 



(0) 



<Mo) 



e„(«) 



K n — 1 



-z6„_i(z) 



Wl(0 le;( z )-^e;_ l( z), 



i(0) 



n* n (z) + n n (z) 



(z6 n ( 2 )-6*(z))+n 



^n-l 

W(.z) 



v+i 



<^n+2(0) , . K n+1 », » 

<Pn+l(0) K n 



W(z) 
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Proof. Consider (|2.66(l first. In Ij2.58(l let us map n ^ n + 1 and add this result to 
(|2.6()|) . However the combination f2„+i(z) + f2* (z), which appears in the resulting 
sum in this form and also under n <— ► n — 1, occurs in 1|2.62[1 allowing it to be 
eliminated. To derive (|2.67() we again map n *—> n + 1 in (|2.58(l and subtract this 
from 1|2.62|) eliminating Sl„ + i(z). If we add (|2.62|1 and (|2.63(1 and simplify the result 
then we arrive at (12.6811. □ 



For a general system of bi-orthogonal polynomials on the unit circle the cou- 
pled recurrence relations and spectral differential relations can be reformulated in 
terms of first order 2x2 matrix equations (or alternatively as second order scalar 
equations). Here we define our matrix variables and derive such matrix relations, 
and this serves as an introduction to a characterisation of the general bi-orthogonal 
polynomial system on the unit circle as the solution to a 2 x 2 matrix Riemann- 
Hilbert problem. 

Firstly we note that the recurrence relations for the associated functions e n (z), 
e* (z) given in (|2.36I2.37|I differ from those of the polynomial systems (12. 2112. 2211 by 

a reversal of the signs of 4> n (0), 4>n(0)- We can compensate for this by constructing 
the 2x2 matrix 

(< I \ e n{z) \ 
<Pn(z) —TT » 

(2.69) Y n (z) := w ) z > 



and note from (|2.39|) that det Y n 



V 

-2z n /w(z). 



>{z)J 



Corollary 2.3. The recurrence relations for a general system of bi-orthogonal poly- 
nomials 1^2.2112.2^) and their associated functions \2.3b\2.3*fy are equivalent to the 
matrix recurrence 



(2.70) 



ri+l 



1 / Kn+lZ 

n-n \</>„+i(0)z 

According to \2.4-l\j the matrix K n has the property det K n = z . 




Corollary 2.4. The system of spectral derivatives for a general system of bi- 
orthogonal polynomials and associated functions $2.52112. 531) are equivalent to the 
matrix differential equation 



(2.71) Y' n 



W(z) 



( 



V 



n n (z) + v(z)-^±±ze n (z) 



Vn+l 



(0) 



fn+1 



(0) 



zQ* n (z) 



n* n (z)~v(z)-^e; l (z) 



Proof. This follows from &i Ml'l SB and employing \±'l\\'l:l'M:Ml.\V?\ . We note 
that TiA n — n/z — w' /w when l|2.67f) is employed. □ 
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Remark 2.3. Compatibility of the relations l|2.7()JI and (|2.71|1 leads to 
(2.72) K' n = A n+1 K n - K n A n , 

and upon examining the (1, l)-component of this we recover the linear recurrence 
(|2.(i3[) . the (1, 2)-component yields (|2.62|) . whilst the (2, l)-component gives (I2.64f) 
and the (2, 2)-component implies (|2.65l) . 

Remark 2.4. There are, in a second-order difference equation such as l|2.36[l or 
(|2.37Jl . other forms of the matrix variables and equations and these alternative 
forms will appear in our subsequent work. Defining 



(2.73) X n (z-t):= 



<Pn+l{z) — — 



<t>n(z) 



w(z) 
e n {z) 



, X*(z;t):-- 



w(z) J 

we find the spectral derivatives to be 



W(Z) 

W > w(z) J 



fa i \ \t( \ -X- W ^ K »^+ 2 (0) . . , 

ir>vA\ urr \v> I M n (z) - V{z) + n —zQ n+l {z) 

(2.74) W(z)X n = z K n+ i4> n+ i(0) 



e n [z) 



-n n (z) - v(z) 



x n , 



(2.75) 

W(z)X% = 



K (z)-V(z) + (n + l)^ ^rZ ^^Ax* 

Z K n+1 (f> n+1 {0) I X n . 



Another system is based upon the definition 



(2.76) Z n (z;t) 



<Pn+l{z) 1-— 



V 



w(z) 

K(z) -% z) 



z n(z;t) 



K(z) - V(z) 



W(Z) 



w(z) J 

and in this case the spectral derivatives are 
(2.77) W{z)Z' n = 

(-w n {z) - v(z) + -^-e;(*) + (n + 1)^ 

_K n+ l Z 

gn+l(0), 



4>n(z) 



w(z) J 



®n+l(z) \ 



\ 



tn+1 



n* n ( z )-v(z)-^e* n (z), 

Kn+l / 
x Z n , 



(2.78) W(z)Z: 
( 



n n (z)-V(z)-^-zG n (z) 

0n+l(O), 



V 



K n +1 



-®n(z) 



K n <f> n+2 {0) Jot f ^ \ 



-z 2 e* n+1 (z) 



-n n (z)-v(z) + 



Kn+l 



-Z® n (z) 



J 
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We end this section with a characterisation of a general system of bi-orthogonal 
polynomials on the unit circle (and their associated functions) as a solution to a 
particular Riemann-Hilbert problem. 

Proposition 2.8. Consider the following Riemann-Hilbert problem for a 2 x 2 
matrix function Y : C — > SL(2,C) defined in the following statements 

(1) Y(z) is analytic in {z : \z\ > 1} U {z : \z\ < 1}. 

(2) on z G £ where £ is the oriented unit circle in a counter-clockwise sense 
and +(— ) denote the left (right) -hand side or interior (exterior) 

(2.79) y + (z) = y_( z ) 

(3) as z — > oo 

(2.80) F(z) = (l + 0(z- 1 )) 

(4) as z -» 

(2.81) Y{z) = {l + 0{z)) 





0(1) OOz™" 1 ) 
0(1) 0(z") , 



is assumed that the weight function w(z) satisfies the restrictions given at the 
beginning of this section. Then the unique solution to this Riemann-Hilbert problem 
is 



(2.82) Y(z) = 



/ 4> n (z) e n (z) 



n > 1. 



Proof. We firstly note from the jump condition (|2.79|) that Yn, Y21 are entire z G C. 
From the (1, l)-entry of the asymptotic condition l|2.80|l it is clear that Yn = Tr n (z) 
is a polynomial of degree at most n. Similarly Y 2 \ = o- n (z) from an observation of 
the (2, l)-component. From the (1, 2)- and (2, 2)-components of the jump condition 
we deduce 

(2.83) Y + i2 — Y-12 = — —iii, Y + 22 — Y-22 = — ^^21, 

z z 

and therefore 

(2.84) Y 12 =(^ W ^\ F22= ^C-(CK(C) 



If 2-7ri£ C, — z J j 27ri£ ( — z 

Consider the large z expansion of Y12 implied by the first of these formulae 

-1 f d C , -2, 



(2.85) Y 12 = -z~ l / _4-«,(C)7r n (0 + 0(z- 2 ). 

Jt 27r < 
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According to H2.80|) the integral vanishes and so ir n (£) is orthogonal to the monomial 
C°. Now take the small z expansion 

(2.86) Yi 2 = XV / ^^CKtCK^ + O^" 1 ). 

From the (1, 2) -component of the condition (|2.81() we observe that all terms in the 
sum vanish and we conclude the 7r n (£) is orthogonal to the monomials C, ■ • ■ , £ n_1 
and the first term which survives has the monomial £ n . Thus ir n (Q ^ <^n(C)i 
and from the explicit coefficient in the (1, 1) -entry of l|2.80|l 7r n (£) is the monic 
bi-orthogonal polynomial <fr n (0/ K n- We turn our attention to F22 and examine the 
small z expansion 

(2.87) 1*2 = / ^-M0a n {C)~ l +O{z n ). 

1=0 J 1 s 

The (2, 2)-component of l|2.8f \ tells us that all terms in the sum vanish and con- 
sequently <7 n (C) is orthogonal to all monomials Therefore <r n (() oc <^>* (C) 
and we can determine the proportionality constant from the (2, 2)-component of 
the asymptotic formula H2.8()(l . by comparing it with 

(2.88) Y 22 = -z- 1 f -^- w (Qa n (C) + 0(z-% 



to conclude <r n (Q — (£). Finally we note that 

(2.89) [ = ±e n (z), fl« = -ie; ( ,), 

when n > 0. We also point out detY(z) = —z n ~ x . □ 



Remark 2.5. Our original matrix solution Y n specified by (|2.69|l is related to the 

solution of the above Riemann-Hilbcrt problem by 

(2.90) 



\ / 1 








Y n (z) = I Q j_ r(z) 2z_ I , y(z) = I k„ I r„(z) 

K„/ \ Ul{z)j 

Our formulation of the Riemann-Hilbert problem differs from those given in studies 
concerning orthogonal polynomial systems on the unit circle with more specialised 
weights |H],0,[n]! [5] and [5]. We have chosen this formulation as it is closest to 
that occurring for orthogonal polynomial systems of the line ^S], the jump matrix 
is independent of the index n which only appears in the asymptotic condition and 
it is simply related to our matrix formulation 



3. Regular Semi-classical Weights and Isomonodromic Deformations 

All of the above results apply for a general class of weights on the unit circle 
but now we want to consider an additional restriction, namely the special structure 
of regular or generic semi-classical weights. 
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Definition 3.1. [|37j] The log-derivative of a regular or generic semi-classical 
weight function w(z) is rational in z with 

(3.1) W{z)w'(z) = 2V{z)w(z), 

where V(z),W(z) are polynomials with the following properties 

(1) deg(W0 > 2, 

(2) deg(V) < deg(W), 

(3) the m zeros of W(z), {z\, z%, . . . , z m } are distinct, 

(4) the residues p k = 2V(z k )/W'(z k ) <£ Z> . 

The terminology regular refers to the connection of this definition with systems 
of linear second order differential equations in the complex plane which possess 
only isolated regular singularities, and we will see the appearance of these later. 
An explicit example of such a weight is that of the form 

m 

(3.2) w(z) = Y[(z- Zj y>. 
Consistent with the requirements of Definition 3.1, we set 

(3.3, fg=E^"- 

j=l j=l J 

To reduce the computational labour and render the ensuing formulae simpler in 
appearance we are going to assume henceforth that one of the finite singular points 
is fixed at the origin, i.e. that W(0) = 0. 

It follows from these definitions that the Caratheodory function satisfies an 
inhomogeneous form of H3.ll) . 

Lemma 3.1 (PO'GH)- Let the weight w(z) be such that w(e 7rl ) = w(e~ vl ). The 
Caratheodory function \2.S(]\) satisfies the first order linear ordinary differential 
equation 

(3.4) W(z)F'(z) = 2V{z)F{z) + U(z), 
where U(z) is a polynomial in z. 

Proof. Following |37| we can write 

(3.5) w(z)F(z) = j + V (CMC) + / ^^[^(*)-^(0MC). 

The last term is a polynomial in z of bounded degree which will be denoted by 
iT3(z) as we are not going to evaluate it explicitly. Differentiating this yields 

(3.6) ±W( Z )F(z) = J «--^(0«,(C) + **(*), 
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and we rewrite the first term of the right-hand side as 



2TTiC((-z) 2 0/ 2lliz J ((-Z) 



2mz J <i£ ( — z 

(3.7) . JL / «g±±m<MO). 



Consequently we have 



' 1 z (lW'(Q + 2V(Q]w(0+w2(z), 



z J 2niQ C, — z 
1 f d( ( + z 



z J 2tti( £ — z 
1 f d( ( + z 



z[W'(z) + 2V(z)]w(() 

{C[W'(0 + m0\ - z[W'(z) + 2V(z)]}w(0 



z J 2iri( i Q — z 
+ n 2 (z), 

(3.8) = [W'(z) + 2V(z)] [ n \i + Z w(C) + ttiW, 



2ni( ( - z 

and (|3.4[) follows. □ 



This lemma leads to the following important set of formulae. 

Proposition 3.1. Assume that W(0) = 0. The coefficient functions Q n (z), 0*(z) 7 
Q n (z), Qn(z) are polynomials in z of degree m — 2, m — 2, m — 1, m — 1 respectively. 
Specifically these have leading and trailing expansions of the form 



(3.9) 9 n (z) = (n + l + V ft )— 



n z m-2 



3 = 1 



K n +1 



S lib lib III, III 



ft:,, 



nl 0„ +2 (O) 
K 2 n+1 K n+2 0n+l(O) 



(0)4(0) K n'n+l I m _3 



Ni+l 



+ 0(z m - 4 ), 
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(3.10) 6„(z) = [2V(0) - nW'(0)] 



0n+l(O) 

«n^n-l(0) 



+ {[2V'(0) - i„W'(0)]-^L + [2 V(0) - (n - 1)W'(0)] 

L 0n+l(O) K„_10 n+ l(O) 

+ f[(n + 1)W(0) - 2V(0)]-^- - [(n - 1)W'(0) - 2V(0)}^-) ^" (0) 



(0). 

+ 0(z 2 ), 



(3.11) e:(z) = -(n + ^ Pj )^-z™- 2 



^ " >n+l(0) 

+ {[(»+ E Pi) f> - E + (» + 1 + E p^^k 1 -^ 

l ~i ~i ~i ^n+ll ) ^ 0„+i(O)K„+l 

n—1 I 

-< n " 1+ 5>> ) 



± [ z m-3 



+ 0(z m - 4 ), 



(3.12) Q* n (z) = -[2V(0)-(n+l)W'(0)} ' 



&n+ 1 

^71 



+ ( - [2V'(0) - |(n + 1)W"(0)]— - [21/(0) - nW'(O)] 
+ [(„ + 2)W(0) - 2V(0)] ( 4 2 ^" +2 ^ - ^In±L] \. 

\«n+2K n+ i <?n+l(0) Kn+lKn+l/J 



m 

(3.13) f7„(z) = (l + -^p,)z m - 1 

m m ^ rn m 

- 2 (E ft)(E ^ + 2 E - E z j 

+ (» + 2 + EPi)— ^— T^-— } 

^ K„ +2 K n+ l n+ l(O) K n+ 1 J 



+ 0(z m - 3 ), 

(3.14) n n (z) = V(0)-nW'(0) 

+ / v"(0) - *nW"(0) + (V(0)-^- + [V(0) - nW'(O)]- 



«n+l\ 0«(O) 



Kn / 0n+l(O) 



+ [V(0) - nW'(O)]— - [V(0) - (n+ l)W'(0)]-^-}« 



+ 0(z 2 ), 
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(3.i5) n;co = 4x>* m -' 



2 ^pf 

3=1 

Kn 4>n(0) , L+l \ m _ 2 



2 h h 2 h ^V+l<^ + i(0) /Sn+J 



2 



+ 0(z m - 3 ), 

(3.16) n*(z) = (n + l)W(0)-V(0) 

+ (i(n+l)W"(0)-F'(0)+[(n+2)W(0)-2y(0)] ^ j" +2 ^ -^'(O)-^L 

I K n+ 2K„+1 0„+l((J) K„ + l J 

+ 0(, 2 ). 

Proof. Following the approach of Laguerre pQ we write F(z) in terms of 4> n {z), 
ipn(z) , e n (z) and use (|3.4|l to deduce 

(3.17) = WF' - 2VF - U, 
= W ( en ~^ n \ _ 2V en -U 



The numerator of the first term is independent of e„, and so is a polynomial in z, 
and we denote this by 

(3.18) 2 ^"+ l(0) z n e n (z) = + e n <#J + 2^„e„. 

Given that this is a polynomial we can determine its degree and minimum power 
of z by utilising the expansions of n ,e n both inside and outside the unit circle, 
namely (|2.4til2.48[l . We find the degree of the right-hand side is n + m — 2 so that 
@ n (z) is a polynomial of degree m — 2. Developing the expansions further we arrive 
at 13.9fl . An identical argument applies to the other combination 

(3.19) 2 ^z»+ 1 e;(z) = w(<&<% - e ;<) - 2v>; e ;, 

and G>n(z) is also a polynomial of degree m — 2 with the expansion i|3.10|l . To 
establish (|3.13(l we utilise the other form of 0„(z) and l|2.38|l to deduce 

(3.20) WM n - <M' n ) + 2V<t>ni>n - U<f n = 2 " +l(O) z"9»(z), 

= [<f> n +\4>n - ll>n+\4>n]®n{z)- 

Separating those terms with <p n and i/>„ as factors we have 

(3.21) {e„(z)0„+i - w<j>' n - v<f> n }i>n = {e„(z)^„+i - wa + v^ n - u<f> n }cj> n , 
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so that this polynomial contains both 4> n and ip n as factors and can be written as 
^n'pntpn with Cl n (z) a polynomial of bounded degree. This latter polynomial can 
be defined as 

2 ^" +l(0) Z n n n (z) = W(i> n+1 <P' n - K+l^'n) + V(M n+1 + Vr^n+l) " U<j> n cj> n+1 , 

(3.22) = W(e n+1 <j>' n - c/)n+ie' n ) + V(c/) n e n+1 + e n <p n+1 ). 



Again employing the expansions 12.46I2.48|) we determine the degree of Vt n (z) to 
be to — 1 and the expansion i|3.13[l follows. Starting with the alternative definition 
of e*(z) and (|2~35|) 

(3.23) w^ir: - «) - 2vr n r n - = 2 ^lz"+ i e;(z), 

and using the above argument we identify for the polynomial O* (z) 

r>>n 

(3-24) = W{-e* n+1 c% + K+i<) " V(K<+i + 

The degree of is found to be to — 1 and it has the expansion l|3.15|l . □ 

Remark 3.1. Solving for 4>' n and e' n between Ij3.18l) and (|3.22(1 leads to (|2.52|l and 
(|2~51> . whilst solving for and e*' using (|3.1i)l3.M|) yields (|2~53|) and (|2~551) . 

Bilinear residue formulae relating products of a polynomial and an associated 
function evaluated at a finite singular point will arise in the theory of the defor- 
mation derivatives later. These are consequences of the workings of the proof of 
Proposition 13. II so we give a complete list presently. 

Corollary 3.1. Bilinear residues are related to the coefficient function residues in 
the following equations, valid for all Zj 

(3.25) Mzj)en{z . )=2 ^±M z - 6 ^ 



(3.26) 4n{*i)<{*i) = - 2 



2v( Zj y 
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(3.27) 

(3.28) 
(3.29) 
(3.30) 
(3.31) 
(3.32) 
(3.33) 
(3.34) 



, , w x O ^n+i(o) n n n ( Zj ) + v{ Zj ) 
^ r, v r„ ^ ^+i(Q) ^ Mz j )-v( Zj ) 

= 2_ ^*' 2^ ' 

o 0n+l(O) „ n+ l ";(^)-n^) 
J 



2T/(z,) 

finfe) - V(Zj) - -^ZjQniZj) 

n* n ( Z3 )~v( Zj )~^e* n ( Zj ] 



il n ( Zj ) + V(zj) - -^-ZjQ n {zj) 



n* n { Zj ) + V{ Zj )-^Q* n { Zj ) 



Proof. These are all found by evaluating one of l|3.18|l . I|3.19|l . (|3.22ll . or (|3.24() at 
z = Zj and using (|2\38l2\40|) . □ 

The initial members of the sequences of coefficient functions {On}.^L i {@n}5£Lo> 
{n n }%L , {K}™ =0 are given by 



(3.35) 
(3.36) 

(3.37) 
(3.38) 



1(0), 
6a(0) 



K 



«1 



%' 2 

n ^2 
3 



z9i(z) = -iz 2 (2F(z) - Kq[/(z)) ~ 2«i0i(O)zE7(» 



,«10l(O) 



W(z) 



#(0) 



(2V(z) + ^CT(z)), 



2 MPi 2 2 e * (z) 
«1 



)=-2V(*)-K§tf(z), 



>i(0) 



z 2 (2^(z) - KoC/(z)) - 2kiM°) zU ( z ) 

_ 2 ^p. w{2) _4 {2V{2) + K 2 U{2))) 



(3.39) 
(3.40) 



20i(O)no(z) =K lZ (2V( Z ) - k 2 U(z)) ~ k 2 q M0)U(z), 
20 1 (O) z n*( z ) = - ki(2V(z) + k%U(z)) - kIM0)zU(z). 



One can take combinations of the above functional-difference equations and 
construct exact differences when z is evaluated at the finite singular points of the 
weight, i.e. W( Z ) = 0. The integration of the system is given in the following 
proposition. 
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Proposition 3.2. At all the finite singular points Zj,j — 1, . . . m, with the excep- 
tion of Zj = 0, the coefficient functions satisfy the bilinear identities 

«n0n+2(O) 



(3.41) 




(3.42) 




(3.43) 





Kn+10n+l(O) 



Zj@ n {zj)@ n+1 {zj) + V (zj), 



"n-1 <Pn+l 



(0) 



i(0) 



0n+l(O)^„(O) 



e„(z,)e;_ 1 (z J ) + ^ 2 (z,), 



(0) 



<Mo) 



ZjQ n (zj), 



(3.44) 
(3.45) 
(3.46) 



(0)^,(0) 



0|e*(^-)e„_ 1 ( 2j ) + ^ 2 (^), 



n+ i(O)^„+i(O) 



ZjQ^e^ + V^Zj) 



.Krsi Proof. We take the first pair of identities Ij3.41|l and H3.42JI as an example for 
our first proof. Multiplying the f2 n , 0„_i terms of I|2.fi8|) by the corresponding terms 
of (|2.59|) , evaluated at a finite singular point z = Zj , one has an exact difference 



(3.47) fln-i(*i) 

K„0 n+2 (O) 



«n-l<^n+l(0) 



z J e n _i(z :; -)e„(z i ), 



K n+ 10 n+ l(O) Kn^n(0) 

assuming none of the Zj coincide with —r n +x/r n for any n. Upon summing this 
relation the summation constant is calculated to be 

K O 02(O) 



(3.48) 



«i0i(O) 

by using the initial members of the coefficient function sequences in 

The result is <|3.41[1 . whilst the second relation follows from an identical argument 

applied to B3SBZBP - □ 

Second Proof. The three pairs of formulae (|3.41I3.42I) . (|3.43I3.44I) and (|3.45I3.46|) 

arise from the fact that at a finite singular point Zj the determinant of the matrix 
spectral derivative must vanish. Thus (|3.41|) and l|3.42|l express the condition that 
the determinant of the matrix on the right-hand sides of l|2.74ll and l|2.75|l vanish 
respectively. It can be shown that the same condition applied to the right-hand 
sides of (|2.78(l and (|2.77(l implies (|3.43() and (|3.44() respectively when one takes 
into account the identities (|2.fi3l) . (|2.(ifi|) . (|2.fi7l) and l|2.fiUf) . The last pair are a 
consequence of det(WA n (zj; £)) — along with the identity I|2.fi7ll . □ 
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Third Proof. All the bilinear identities in Proposition ^ . 2l can be easily derived from 
the residue formulae 1(3.2513.34(1 by multiplying any two of the above formulae and 
then factoring the resulting product in a different way. Thus 1(3.41(1 arises from 
multiplying 1(3.27(1 and 1(3.28(1 and then factoring the product in order to employ 
((3.25(1 . Equation 1(3.43(1 comes from multiplying 1)3.25(1 and 1(3.26(1 with nnn-1, 
using the recurrences 1(2.22(1 . 1(2.37(1 with n i— > n — 1 to solve for 0*„ 1 (z J ), tn-ii z j) 
and employing (|3.25l) along with l|3.27|l and H3.28|l setting n i— > n — 1. Equation 
(I3.45|l is derived by multiplying 1(3.25(1 and 1(3.26(1 and then factoring using (|3.31() 
and 13.33fl . The reciprocal versions follow from similar reasoning. □ 

Remark 3.2. It is clear from the first proof that the bilinear identities given in 
Proposition 13.21 can be straightforwardly generalised to ones that are functions of 
z rather than evaluated at special z values. They can be derived directly from 
Proposition 12.71 so apply in situations where the weights are not semi-classical, 
and contain additional terms with a factor of W{z) and sums of products of other 
coefficients ranging from j = 1, . . . n. However because we will have no use for such 
relations we refrain from writing these down. 

Remark 3.3. As z = is a finite singular point then the limit as z — ■> may be taken 
in the product of ((2.5812.59(1 . however this does not lead to any new independent 
relation but simply recovers 



In the case of a regular semi-classical weight function the matrix A n (z;t) has 
the partial fraction decomposition 



under the assumptions following ((3.1|) . Let us take the first finite singularity to be 
situated at the origin, z\ = 0. Some care needs to exercised as many relations differ 
depending on whether Zj = or not because of the additional term in ((2.67|1 . The 
residue matrices for the finite singularities are given by 



Q n (0) = 7(0) - nW'(0). 



m 



(3.49) 




(3.50) A 



nj — 



2V(z 3 ) 



X 



/ 



V 



n n { Zj )-v{z j ) + 

0n+l(O) 




for j — 2, . . . , m and 



(3.51) 



Anl 



27(0) 



Pi 



nW'(O) - 27(0) [27(0) 




nf'(0)]r, 
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Using the identity 1)2. 67|) we note that 

(3.52) TrA n j = -p j} j = 2, . . . ,m, 

(3.53) TrA nl = n- p x . 

In either case we find that det A n j = using l|3.45[l . An alternative expression for 
the residue matrices in the case j = 2, . . . , m is 

(3 54) An- = --Pjz- n ( ^ z ^ €n ^ -Mzj)tn(z 3 )\ 

2 3 3 \-<l>n(Zj)£n(Zj) <t>n{Zj)z n { z i) ) ' 

The regular singularity at z — oo has a residue matrix given by 

m 

(3.55) A noc :=Res z = {-z^Aniz- 1 )} = -^A nj . 

Using this definition and the large z terms for the coefficient functions l|3.9|l , l|3.11|l , 
(|3.I3|I . I|3.15(l we evaluate this matrix to be 

(3-56) A noc = [ ] . 

We read off that TrA noo = — ti+J^ - =1 /Oj and det j4 noo = — n X^j=i Pi- Furthermore 
relations 1)3. 56[) and (|3.55l) imply the summation identities 



^ m 

(3.57) - ^ PjZj n <i>n{zj)en{zj) = 0, 

i=i 

^ m 

(3.58) - PjZj n <t>n{ z j) e n{zi) = ~n, 
(3-59) - P] zj n Mz )el{z ) = £ Pi , 

1 m m 

(3-60) - £ ft^CfeKfe) = (n + £ ft)*V 

i=i j=i 

We wish to close this part by commenting on how one would obtain the discrete 

analogs of the Schlesinger equations or multi-variable extensions of the discrete 

Painleve equations from the theory outlined above. If one evaluates (|3.41() (or 

(|3.42|l for that matter) at two distinct singularities 21,2:2, consolidates terms and 

then takes their ratio the result is 

zie„(zi)e n+ i(zi) = [» n (z 1 )-u(z 1 )][r>»(z 1 ) + u(z 1 )] 

Z2Qn(z 2 )Qn+l(z 2 ) [fi„ (*&) ~ V(z 2 )} [Sl n {z 2 ) + V(z 2 )} ' 

This constitutes a recurrence relation for Q n (zx)/Q n (z 2 ). To find a recurrence 
involving f2„ one adopts another method. By comparing the expansions inside 
and outside the unit circle, for example equations l|3.9|l and l|3.10|l for 0„(z), and 
in particular where they overlap one can derive expressions for the sub-leading 
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coefficients l n ,m n in terms of the higher ones. Employing these expressions in one 
of the expansion forms for 0^(^1,2) f i|3.13|) or l|3.14[l ) will yield a recurrence for f2„. 

We now consider the dynamics of deforming the semi-classical weight (|3.3|) 
through a ^-dependence of the finite singular points Zj(t), 

m 

(3.61) » = _^ Pi _S_ > 

W ^ Z — Zj 

3=1 J 

where ' := d/dt. Given this motion of the finite singularities we consider the 
i-derivatives of the bi-orthogonal polynomial system. 

Proposition 3.3. The deformation derivative of a semi-classical bi-orthogonal 
polynomial is 



m . m 

(3.62) Mz) = {-^-J2pA + - Y.pA^niz.m^^Uniz) 

^* t\ijy Zn £s Z 1 2s 2/ ' n * 

j — 1 J j — 1 J J 

m -. 

v £i Z'i Z Z j ■» 

3=1 J J 

whilst that of its reciprocal polynomial is 

(3.63) fa(z) = { - ^ + ^Eft?^ n <fe)0nfe)-— 

/'VtT — 2s 7 2s 2s 7 s> 

3 =l J J 

^ rn 

~ { 9 Pi-t. Z T e n i^Wn ( Z 3 ) " T I ^" • 

3 = 1 3 3 

The deformation derivative of an associated function is 

-. rn 

(3-64) e n (z) = { - -£ Pi 7^- }*»(*) 

3 = 1 J J 

1 " l 1 

+ loE' , c z r" 6 ''( z ^ ,i fe); — 

v di Z'i Z Z j * 

3 = 1 

and that of a reciprocal associated function is 



-. rn 

(3.65) e* n (z) = { - ^ - -^^zj-^felCfe) — 

+ { 9 S ft S Z 7" e « ^ ^ ) T~~T } £ « ( Z ) • 

^ 2-7" Z Z j J 

3=1 

Proof. Differentiating the orthonormality condition 



2ni( 
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and using (|3.61(l we find 

0=— ■z^w(04>n<f>n-i-y~]pjZj [ ■^rrw(C)—^ — 0„0„_i, i = 0, ..../;. 

K n J 2tti( ; J 2m C C - Zj 



Now 



:w{Q- <p n (0<f) n -i(0 = / — — zw(Q(p n (Q- 



. 7 / -In f d C ,^ M0 

= -h,o+ * n - i } Z ~ 1) e n (z), n>0, 
z 2z 



so that 



In addition we can represent <p n -i(z) as 



-t0) =^2^i,j ( t > n-j{z), 



3=0 



3=0 



Yl I -^^ W (0 ( t>n-i{0 ( t>n-j(0 ( t>n-j(z), 



n _ 7 ' ), 



S 1=0 

:w(Q(t>n-i{()- 



2iri( WT w l-(z 
Writing the Kronecker delta in a similar way the whole expression becomes 



-3=1 ' '-^i' 

for all < i < n and H3.(j2|l then follows. The second relation follows by an identical 
argument applied to 

The derivatives of the associated functions l|3.64|) , l|3.65[) follow from differentiating 
the definitions 1)2. 32|) , (|2.33|l and employing the first two results of the proposition 
along with the relation (|2.40l) . □ 
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Corollary 3.2. The t- derivatives of the reflection coefficients are 

Proof. An alternative formula to i|3.62[l is 
(3.68) 

(m . \ 1 m n 

and by examining the coefficients of ( n , (° we deduce that 

1 m 

Noting that the derivative term of (|3.22|) vanishes when z = Zj and employing i 
we arrive at l|3.66|l . The second equation, (|3.b7[) . follows by identical reasoning. □ 

Sums of the bilinear residues over the finite singular points are related to defor- 
mation derivatives in the following way, 

rn . m 

2 , = _^ Pj z -± + ^E^r^^te), 

.7=1 J .7=1 J 

^ m 

(3-69) = - »^Pi^-zj n e* n (zj)4> n (zj), 



■ ,_, m . m 

0n(Oj K n V - ^ z j _ 1 K n ST^ z 



E^f = oZTnrE^f z j 



<M«) K„ ^ *j 2^„(0) ^ 



0„(O) k„ Ik, 



EpiS*r<(*i)&&)' 



</>„(o) k„ 2 0„(o) fri zj 

(3-71) = ^# V -^- Zj QUz 3 ). 

V ' ^(0) ^ 2V( Zj ) 3 nK 31 

For the regular semi-classical weights we can also formulate the system of de- 
formation derivatives as a 2 x 2 matrix differential equation and demonstrate that 
the system preserves the monodromy data with respect to the motion of the finite 
singularities Zj(t). 
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Corollary 3.3. The deformation derivatives for a system of regular semi- classical 
bi-orthogonal polynomials and associated functions IS.fiOMS.fity are equivalent to the 
matrix differential equation 



(3.72) 



Boo ^2 _ _V ^nj \ Yn 



3=1 



whe 



(3.73) 



( *± \ 



\ 



Proof. This follows from a partial fraction decomposition of the system I|3.62I3.65|) 
and using (|3.6i)l3.7lj) . □ 

Remark 3.4. For a special class of irregular semi-classical weights a result analogous 
to Corollary 13. 31 has been given by Bertola, Eynard and Harnad |15j . 

In the case of the pair l|2.70|l , (|3.72|) compatibility implies the relation 



(3.74) 



K n = B n+ iK n - K n B n , 



however there are no new identities arising from this condition. Taking the (1, 1)- 
component of both sides of this equation we see that it is identically satisfied 
through the use of i|2.63[) and i|3.71() . Or if we take the (1, 2) -components then 
they are equal when use of made of l|2.58|l and I|3.69I3.66[1 . In a similar way we find 
both sides of the (2, l)-components are identical when we employ (|2.tj(J|) and l|3.71ll . 
Finally the (2, 2)-components on both sides are the same after taking into account 
gB3 and <3.(i9l3.7lt . 

For the pair of linear differential relations l|2.71|l . I|3.72|l compatibility leads us 
to the Schlesinger equations 



(3.75) 
(3.76) 



Anj — [Boo, A n j] + S ' 



Zk 



Zj Zk 



Again there is not anything essentially new here, that couldn't be derived from the 
system of deformation derivatives I|3.62I3.65J) , but it is an efficient way to compute 
the deformation derivatives of bilinear products. Employing the explicit repre- 
sentations of our matrices A n j we find the following independent derivatives in 
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component form 

d Pj 



(3.77) 



dt2V(zj) 



n n ( Zj ) + v(zj)- 



"n+1 



ZjO n (zj) 



_ Pj 0n+l(O) d 

2V(zj) k% dt 
Pj 0„+i(O)0„+i(O) ^ p k Zj - Zk 

2V(z k ) zj - z k 



2V( Zj ) 

(3.78) 
d pj 



(K n 4> n (0))Q n (zj) 

[z k @l{Zk)®n{Zj) - ZjG n (z k )@n( z j)] 



(0) 



dt 2V(zj) K n 



<dn{Zj 



Pj 4>n+l(0 ) / «• 



-e„(z,) + E 



21/(z fe ) Zj - z fc 



Q n (z k )[£l n (zj) 



Kn+l 



ZjO„(zj)] - Q n (zj)[Cl n (z k ) 



Z k Qn(zk)] 



(3.79) 



Pi 



,J n+l 



(0) 



2j n(*j) 



'"' +l(0) J - ?*e;(*,) 



:(«n^(0))[n;(^) 



Kn+1 



<*;&)]-£■ 



« fc e;(^)[n;(«i) - — e;^)] - ^(^[n^) - ^e;(z*)] 

The fact that the deformation equations satisfy the Schlesinger system of partial 
differential equations should be of no great surprise as the isomonodromic properties 
of the regular semi-classical weights are quite transparent. 

Proposition 3.4. The monodromy matrix Mj, j = 1, . . . , to, oo defined by the 
analytic continuation of Y n around a closed loop enclosing the singularity Z j 

(3-80) Y n \ z . +5eM = Y n \ Zj+g Mj, 

is constant with respect to the deformation variable, Mj = 0. 

Proof. In the neighbourhood of any isolated finite singularity \z — Z j\ < A the 
Caratheodory function is a solution of the inhomogeneous first order differential 
equation 1)3. 4)1 and can be decomposed as 



(3.81) 



F(z) = fj(z) + Cjw(z). 



Here fj(z) is the unique, holomorphic solution of the inhomogeneous ODE in this 
neighbourhood whose existence is guaranteed by conditions (3) and (4) of the def- 
inition 1)3. ljl and which we express as 



(3.82) 



l>0 
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The second term of 1(3. 81(1 is a solution of the homogeneous form of 1(3.4(1 and Cj a 
coefficient depending only on {zk, Pk}kLi- We note that the deformation derivative 
of the Caratheodory function can be written as 

(3.83) Hz) = ^f(z) + \Y,p^ z ^n^) - 5«*5>J. 

under the assumptions that either z^ ^ T, k = 1,..., to or if any z^ e T then 
SRpfe > 0, to ensure the integral l|2.3Q(l is uniformly and absolutely convergent. 
Differentiating (|3.81(l with respect to t then implies the relation 



(3.84) Cj(z - ztf* f[(z - z k r + Pj -%— \f S {z) - F {zi)\ 

~ y 3 



+^fe-(/+l)^/+i)(z-^)'+^p/— ^[/,(z)-^(z/)]+^p^K-F(z ; )] = 0. 

l>0 tyj I 

Making the assumption dtpj > we can use the equality F(zj) = fj(zj) and 
consequently under condition (3) of 1)3.1(1 we have 

(3.85) Cj( z ~ z j) Pj + analytic function of z = 0, 

for \z — zj\ < min{A, \zj — Zf.\, \\ z j\ — This implies Cj = by condition (4). 
Using this fact and that the monodromy matrix is given by 

<w "> - (S Ci(1 e : s C") • 

the proof is concluded. □ 

Remark 3.5. The monodromy matrices are all upper triangular, which is consistent 
with the fact that we are dealing with classical solutions of the Schlesinger systems. 
Also, they are independent of n, which is to say they are preserved under the 
iteration n — ► n + 1. 
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